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Abstract 

The equivalence is established between the one-dimensional (ID) 
Bose-system with a finite number of particles and the system obeying 
the fractional (intermediate) Gentile statistics, in which the maximum 
occupation of single-particle energy levels is limited. The system of ID 
harmonic oscillators is considered providing the model of harmonically 
trapped Bose-gas. The results are generalized for the system with 
power energy spectrum. 
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1 Introduction 

The observation of Bose-Einstein condensation (BEC) in ultracold trapped 
alkali gases [H [2] gave new stimulus to the study of this quantum phe- 
nomenon. In particular, the effects of a finite number of particles on BEC 
of an ideal gas was discussed in theoretical works [31 HI El El [71 [8] , where the 
corrections to the bulk properties were found. 
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The aim of the present study is to propose the description of a bosonic 
system with a finite number of particles by means of finding such a model 
system, for which the treatment might be mathematically simpler. To some 
extent, such an approach has common features with finding boson-fermion 
equivalence in ideal gases [9] or Tonks-Girardeau gas [10] achieved experi- 
mentally in 2004 The anyon-fermion mapping is also known in the 
application to ultracold gases [13]. 

Haldane's exclusion statistics [H] was considered by Bergere [15]. The 
connection of the exclusion (anyon) statistics parameter and the interaction 
in one-dimensional systems was studied in [16], [TTl [18] . recently, a combina- 
torial interpretation of exclusion statistics was given by Comtet et al. [T^ . 

Another approach is seen in a different type of the fractional statistics, 
which is formally understood as an intermediate one between Fermi and Bose 
statistics. Namely, the maximum occupation of a particular energy level is 
limited to M, with M = 1 corresponding to the fermionic distribution and 
M = oo being the bosonic one, respectively. This statistics is known as the 
Gentile statistics [201 EB [22] • If the relation between the number of particles 
N in the real system and the parameter M in the model one can be found, 
the stated problem is solved. 

For simplicity, a one-dimensional system is considered. The paper is 
organized as follows. Microcanonical approach for harmonic oscillators with 
single-particle energy levels given by Em = hum is considered in Sec. [21 The 
oscillators, unlike classical particles, are indistinguishable reproducing thus 
a quantum case. Physically, this corresponds to bosons trapped in a highly 
asymmetric harmonic trap. In Sec. [3l the same system is treated within 
canonical and grand-canonical approaches. Sec. [H contains the generalization 
of obtained results for the system with power energy spectrum Em oc m*. 
Short discussion in Sec. \5\ concludes the paper. 

2 Microcanonical approach 

The number of microstates r{E) in the system of ID oscillators is the num- 
ber of ways to distribute the energy E = hujn over the (indistinguishable) 
particles. Such a problem reduces to the problem in number theory known 
as the partition of an integer [231 [211 [25]. An asymptotic expression for (un- 
restricted) partition is given by the well-known Hardy-Ramanujan formula 
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Thus one obtains: 
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the following equation of state is obtained: 



E = —hu — 
6 \hujj 



As the energy E is an extensive quantity, E (x N, where N is the number 
of particles, the thermodynamic limit uN = const follows immediately from 
the above equation. The same result also might be obtained from different 
considerations [27j . 

If one considers a finite system of bosons or a system of particles obeying 
fractional statistics the number of ways to distribute the energy E = hujn over 
N particles is the problem of restricted partitions of an integer number n [25] . 
For convenience, hereafter huo is the unit of both energy and temperature. 

The expression for the finite system is given by the number of partitions 
of n into at most summands and asymptotically equals [28] : 




The result reducing to the fractional statistics was considered by Srivatsan et 
al. [29], it corresponds to the number of partitions of n where every summand 
appears at most M times: 



Comparing the entropies S'fin = In / fin and Sfrac = In /"frac, one finds the 
equivalence condition linking the maximum occupation parameter M and 
the number of particles A^: 
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3 Canonical and grand-canonical approach 

It is straightforward to show that in the case of the defined fractional statistics 
the occupation number of the energy level e equals [20l [211 [22] 

jM{e,fi, I ) - _ ^ g(M+l)(e-M)/T _ 1' 

where n is the chemical potential and T is the temperature. 

The chemical potential is related to the number of particles A/" as follows: 



^f = J2fMiej,^^i,T) (8) 

j=0 



and energy E equals 



E = Y.e,fM{ej,fi,T). (9) 
i=o 

However, the case of a finite system is much easier to implement in the 
canonical approach. It is possible to show that the partition function of 
indistinguishable ID oscillators is given by (cf. [21]): 

ZN=u{^-e'^^y > (10) 

from which the energy Ej^ can be calculated. In the limit of large the 
leading term is given by 

En — Ebosb ~ gi/T _ (ll'' 

where E-Bose is the energy of an infinite bosonic system. 

For the fractional-statistics system the grand-canonical approach is used. 
The fugacity z = e^/^ is represented as z = zbosg + with zbosg satisfying 

A^-E ,- (12) 

It is found that Az ~ — in the limit of large M, from which the correction 

M 

to the energy given by Eq. ([9]) follows: 

1 

Em — Ebosc ^ J^- (l^) 



Comparing Eqs. fllip and (fT^ one obtains the following relation between the 
parameters M and N: 

M ~ le^/^. (14) 

In the exponent, the temperature T is related to the energy level n of (EI) via 
Eq. ([3]) (with E = n). Result (IT^ thus reproduces the microcanonical one 
([6]) up to the negligible factor of 1/A^ — it must be taken into account that 
only leading terms were preserved in the logarithms of (jl]) and ([5]). 



4 Power energy spectrum 

In this section, a general power energy spectrum e = am'^ (s > 0) is consid- 
ered. By choosing appropriate energy units, one can set the constant a = 1. 
In fact, only s = 1 and s = 2 cases are realized in real physical systems [2H], 
but other values can effectively occur in some exotic model systems. 

To obtain /fin(^) for arbitrary s it is worth to recall briefly the derivation 
of the expression for restricted partitions from [25j . 

Partition function Z{(3) and the number of microstates r[E) are related 
via the Laplace transform: 



/■oo 1 c+ioo 

Z(/3) = / r{E) e-P^ dE, r{E) = — / Z{(3) e^^ d^. (15) 
Jo Zm J-ioo 

The entropy S{P) equals 

S{p) = (3E + In Z{(3). (16) 
For energy spectrum Em = rn^ the partition function is 

m=I[ (l-e-^"T . (17) 



m=l 



Using the saddle-point method, one can evaluate r{E) (ITS!) as follows: 

riE) = ^^EEiM. (18) 
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The entropy S{(3), after applying the Euler-Maclaurin summation for- 
mula, can be expressed in such a form 

SiP)=pE-f:\n{l-e-^^')=PE + ^ + llnP + ..., (19) 
m=i y ^ 

where 

c{s) = r (i + c (i + ^) , (20) 

T{z) and ({z) being Euler's gamma-function and Riemann's zeta-function, 
respectively. 

The stationary point Po is 

f3o = {^^y^^^'^ = XsE~^/(^^'\ (21) 
Thus, the number of microstates is 



r{E) = exp [a.(. + 1)E^ . (22) 

Substituting E with n one can obtain the well-known Hardy-Ramanujan 
formula [26] for the number of partitions of an integer n into the sum of sth 
powers. 

When the number of particles N in the system is finite, the correction to 
the above formula must be found. In this case, the partition function equals 

lnZ;v(/3) = -Eln(l-e"^"'0 (23) 

m=l 

and for the entropy one has 

N 

SUP) =PE-Y.ln{l- e-^"^') . (24) 

m=l 

After simple transformations it is easy to obtain the following: 

SUP) = s{p) - [Kpn^) , (25) 
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where r(a, x) is incomplete r-function. Thus 
rUE) = r{E)exp ^ 



(26) 



Applying the asymptotic expansion for r{a,x) [5U^, Eq. 6.5.32], we finally 
arrive at the following: 



rfi„(^) = riE) exp 



s(3o 



-N'~'e 



(27) 



Substituting E with n one obtains the result for restricted partitions 



rfin(n) = r(n) exp 



-n 



/{s+l)j^l-s^-XsN''n-''/(''+^^ 



(28) 



cf. also Eq. (17) from [31]. For this function, the notation p%{n) is tradi- 
tionally used, note, however, that in the problem of integer partitions s must 
be integer. For s = 1 the obtained expression reduces to that of Erdos and 
Lehner [2B], see Eq. (Hj). 

The fractional-statistics result can be directly taken from 
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To obtain the relation between the parameters M and N , one can again 
consider the entropies S'^ac = InTfrac and S'gn = In /fin: 



5frac -S = -^iVi-e-'3o^\ S^^-S = -^M-'l% (30) 
s(3q 1 + s 

where S = \nr. 

Dropping the constants, the following result is obtained: 

M^/" ~n^iV*-^exp{A,n"^A^'}. (31) 
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It is interesting to find in this general case the connection between energy 

X dS 

E and temperature T from the definition — = : 

^ T dE 

^ = XsE-^s ^ = A,T. (32) 

Thus, the leading contribution in the relation of M and (1311) is 

M-expU- , (33) 



T 



which is compatible with (fT^ . 



5 Discussion 

To summarize, the equivalence is established between the finite bosonic sys- 
tem and the system obeying fractional (intermediate) Gentile statistics in 
the case of one-dimensional harmonic trap. This approach is extended to a 
general power energy spectrum. While the expressions for two-dimensional 
(2D) partitions are also known [23], the application to asymmetric (elliptical) 
traps as well as the generalization for arbitrary 2D systems needs additional 
study. 

Interacting systems are of special interest now. Weak interactions are 
known not to change the properties of a Bose-system drastically. Thus, 
one can use, e. g., a slightly modified excitation spectrum [32j and, upon 
calculating the properties of a model fractional-statistics system, obtain the 
results for a finite one from the established equivalence. 
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